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ABSTRACT 

I have re-visited the spatial distribution of stars and high-mass brown dwarfs in the 
a Orionis cluster (~3Ma, ~360pc). The input was a catalogue of 340 cluster mem- 
bers and candidates at separations less than 30arcmin to a Ori AB. Of them, 70% 
have features of extreme youth. I fitted the normalised cumulative number of objects 
counting from the cluster centre to several power-law, exponential and King radial 
distributions. The cluster seems to have two components: a dense core that extends 
from the centre to r « 20arcmin and a rarified halo at larger separations. The radial 
distribution in the core follows a power-law proportional to r^, which corresponds 
to a volume density proportional to r"^. This is consistent with the collapse of an 
isothermal spherical molecular cloud. The stars more massive than 3.7 Mq concen- 
trate, however, towards the cluster centre, where there is also an apparent deficit 
of very low-mass objects (M < 0.16 Mq). Last, I demonstrated through Monte Carlo 
simulations that the cluster is azimuthally asymmetric, with a filamentary overdensity 
of objects that runs from the cluster centre to the Horsehead Nebula. 

Key words: open clusters and associations: individual: a Orionis - stars: formation 
- stars: low mass, brown dwarfs. 



1 INTRODUCTION 

The (T Orionis region in the Ori OB 1 b association is fi- 
. nally becoming recognised as one of the most important 
' young open clusters, with an age of only about 3 Ma. In 
the discovery paper. Garrison (1967) used the term "clus- 
tering" to refer to an agglomeration of fifteen B-type stars 
surrounding and including the multiple star a Ori. After- 
wards, Lynga (1981) tabulated a Orionis in his catalogue of 
open clusters. Since the rediscovery of the cluster by Wolk 
(1996) and its subsequent study in depth, which has revealed 
the most numerous and best known substellar population 
(Bejar et al. 1999; Zapatero Osorio et al. 2000, 2002; Ca- 
ballero et al. 2007), only a few authors have investigated 
the G Orionis spatial distribution. In particular, Bejar et al. 
(2004) and Sherry et al. (2004) analysed the radial distribu- 
tion of (7 Orionis cluster members and candidates in annuli 
of width Ar as a function of the separation r to a Ori AB. 
To maximise the number of objects per annulus and min- 
imise the Poissonian errors, Ar must be wide. This leads to 
have few annuli (no more than 12 in the r — 0-30 arcmin 
interval) to fit to a suitable radial profile (exponential de- 
cay - Bejar et al. 2004; King - Sherry et al. 2004). Both 
studies agree that the cluster may extend only up to ~25- 
30 arcmin. The low surface density of cluster members at 
larger separations, the sharp increase of extinction due to 



the nearby Horsehead Nebula-Flame Nebula-IC 434 com- 
plex and the closeness to (or even overlapping with) other 
stellar populations in the Orion Belt surrounding Alnitak 
{( Ori) and Alnilam (e Ori) prevent from suitably broaden 
the radial distribution analysis (Caballero 2007a). At the 
canonical heliocentric distance to a Orionis of 360 pc (e.g. 
Brown, de Geus & de Zeeuw 1994), the cluster would have 
an approximate radius of 3 pc. 

In spite of the agreement on the size of a Orionis, the fits 
and the profiles in Bejar et al. (2004) and Sherry et al. (2004) 
seem to be rather incomplete and inappropiate, respectively. 
On the one hand, the King models were designed for tidally 
truncated globular clusters (King 1962, 1966; Meylan 1987), 
and have also been satisfactorily used for describing galaxies 
(e.g. Kormendy 1977; BinggeU, Sandage & Tarenghi 1984). 
These systems have had enough time to be isothermal, on 
the contrary to very young open clusters like a Orionis, 
where only gravitational relaxation by initial mixing may 
have occurred (King 1962). On the other hand, Bejar et al. 
(2004) exclusively focused on the cluster substellar popu- 
lation. Besides, the exponential fit in Bejar et al. (2004) 
only accounted for the five innermost annuli, which leaded 
to a high uncertainty in the derived parameters. Last, in 
the works by Sherry et al. (2004) and Bejar et al. (2004), 
the input list of cluster members and candidates came from 
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Figure 1. Binned surface density profile for cr Orionis witli Ca- 
ballero (2007c) data. The mean number of objects per annulus is 
34 for a widtli Ar = 3 arcmin. All the figures are in colour in the 
online article on Synergy. 

VRI / IZ optical surveys. Many sources in both analysis had 
no near-infrared or spectroscopic follow-up. 

For a correct study of the spatial distribution in a Ori- 
onis, it is therefore necessary to use new fitting radial pro- 
files and an input catalogue as comprehensive as possible. 
It must cover a wide mass interval. Maximum complete- 
ness and minimum contamination of the catalogue are also 
desired. These requirements are verified by the Maynt cat- 
alogue, which tabulates 339 a Orionis members and candi- 
dates in a 30 arcmin-radius circular area centred on a Ori AB 
(Caballero 2007c). Of them, 241 display features of extreme 
youth (e.g. OB spectral types, Li I in absorption. Ha in 
strong emission, spectral signatures of low gravity, near- and 
mid-infrared excesses due to discs). The catalogue covers 
three orders of magnitude in mass, from the ~18+12M0 
of the O9.5V-hB0.5V binary a Ori AB to the -0.033 Mq 
of the brown dwarf BOS 2.03-617 (Caballero & Chabrier, 
in prep.). Accounting for a Ori A and B as different objects 
separated by ~0.25 arcsec, then the equatorial coordinates 
of 340 young stars, brown dwarfs and cluster member can- 
didates are available. I will use this input catalogue to in- 
vestigate the radial and azimuthal distribution of objects in 
the a Orionis region. 



2 ANALYSIS AND RESULTS 

2.1 Cluster centre and radial gradient 

Caballero (2007a) showed that —46 % of the mass in the 
a Orionis stars with M > 1.2 Mq is contained in the quintu- 
ple Trapezium-like system that gives the name to the cluster, 
and —29 % only in the AB components. This suggests using 
the binary as the cluster centre (r = 0). From the masses 
for the objects in the Mayrit catalogue derived in Caballero 
& Chabrier (in prep.), one third of the total cluster mass 
is encircled in the innermost 5 arcmin. If the mass were ho- 
mogenously distributed within the survey area, the inner- 
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Figure 2. Euclidean minimum spanning tree for a Orionis (Q ~ 
0.88). I have used a matlab/octave function based on an open 
code by F. W. J. van der Berg (University of Copenhagen), based 
in its turn on an algorithm by Hillier and Lieberman (2001). The 
minimum spanning trees for p Ophiuchus, IC 2391, IC 348, Taurus 
and Chamaeleon are shown in Cartwright & Whitworth (2004). 

most 5 arcmin would contain only 2.8 % of the total cluster 
mass [(5/30)^ « 0.028]. I will consider a Ori AB as the ori- 
gin of coordinates because of: (i) simplicity (the coordinates 
of the binary are well determined by Hipparcos; the actual 
coordinates of the cluster barycentre may change when a dif- 
ferent input list of cluster members and individual masses 
is used) ; (ii) reflection of the geometry of the Mayrit survey 
in Caballero (2007c), which was centred on a Ori AB; and 
(Hi) uniformity with previous works (especially with Bejar 
et al. 2004 and Sherry et al. 2004, who also used a Ori AB 
as the coordinate origin). There might be an additional rea- 
son: the largest mass aggregation is probably associated to 
the densest region of the original molecular cloud where the 
fragmentation and star formation initiated (assuming that 
the origin of the referece frame is locked to the -moving- 
cluster barycentre). This reason may be unconvincing, be- 
cause highly turbulent (and fractal-like?) molecular clouds 
probably do not have a "centre" that can be deflned in any 
sensible way, as shown in the simulations of Bonnell, Bate 
& Vine (2003). See Section [23] in this work and fig. 1 in Ca- 
ballero (2007a) for pictorical views of the spatial distribution 
of confirmed and candidate cluster members in the a Orionis 
region, and Caballero (2007b) for a description of the clus- 
ter centre. The old-fashioned plot of the surface density is 
shown in Fig. [T] see similar plots in Sherry et al. (2004) and 
Bejar et al. (2004) for comparison. 

There exists a Q-parameter which quickly, and simply, 
shows if a distribution of cluster members is smooth (a large- 
scale radial density gradient) or clumpy a (multiscale frac- 
tal subclustering; Cartwright & Whitworth 2004). The Q- 
parameter is defined by: 



where m and s are the edge length of the Euclidean mini- 
mum spanning tree and the separation between cluster mem- 
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bers, respectively. A minimum spanning tree is a network 
("graph") of A'max-l lines ("edges") that connect the A''niax 
objects ("nodes") in the shortest possible way under the 
condition no closed loops allowed. See, e.g., Graham, Clowes 
& Campusano (1995) for an application of minimum span- 
ning trees in Astrophysic^. The normalization factors are 
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for m and s (A'^„ 



of the circular survey). 

For the 340 a Orionis cluster member and candidates, 
I have measured m = 0.589 and s = 0.668. Therefore, the 
Cartwright & Whitworth (2004) parameter is Q 0.88. 
This value is larger than 0.80, which distinguish g Orio- 
nis as a cluster with a smooth large-scale radial density 
gradient and a moderate degree of central concentration. 
This concentration is larger than in p Ophiuchus, other clus- 
ter with a radial density gradient, but less than in IC 348 
(Cartwright & Whitworth 2004). Other sparse clusters and 
star-forming regions, like IC 2391, Taurus and Chamaeleon, 
have Q-parameters in the interval 0.47-0.67, which indicates 
that they have, on the contrary, subestructure with fractal 
dimensions between 1.5 and 2.5. 



2.2 Surface density and cumulative number of 
objects 

I present an innovative, accurate, simple method to derive 
the actual expression of the surface density as a distance 
from the cluster centre, o"(r). It can be applied to other 
open and globular clusters and galaxies. The normalised cu- 
mulative number of objects counting from the cluster centre, 
/(r), is: 



fir) 



N{r) 

A^(rma: 



(2) 



where N{r) is the total number of stars in projection within 
a distance r of the centre. If there is azimuthal symmetry, 
N{r) is related to the surface density through the following 
expression: 



(3) 



N(r) = 2tv dr r a[r ) 



The relatively high value of the Q-parameter of a Orionis 
supports the hypothesis of azimuthal symmetry in this clus- 
ter in particular. For systems without azimuthal symmetry 
(e.g. elliptical galaxies), use instead: 

N{r)= / / dr'de r' a(r' ,e). (4) 
Jo Jo 

The surface cr(r) and volume p{r) densities are linked 
through the simple relation a{r) = 2rp(r). This equality 
comes from: 

N(r) = 4n I dr' r'^ p{t') = 2n I dr' r' cr{r') [Vr], (5) 
Jo Jo 



^ The first algorithm for finding a minimum spanning tree was 
developed to find an efliicient electrical coverage of Czecli Moravia 
(Boruvka 1926) 



assuming again azimuthal symmetrjQ. The function /(r) 
varies from at r = to 1 at r = rmax. In the Mayrit 
survey, r^ax = 30 arcmin and A^(rmax) (= A''max) = 340. In 
the discrete approximation, N{r) ~ N*(r) and: 
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I have investigated several functional expressions of a(r) 
that fit in more or less degree the observed normalised cu- 
mulative number of objects, f*(r). A general expression for 
a power-law surface density of index 5 — 2 is: 

5N 

'^(^''^) = ^ri3 ' (7) 

■^n I max 

which, after integration, leads to a simple expression 
for /(r): 

,5 

(8) 



/(r,5) = (^) 

V ^max / 



In this approach, the objects are uniformly distributed in a 
circular area if 5 = 2 (a = constant). Surface densities with 
parameter 5 < 0, that predict a lower number of objects 
close to the centre, were obvioulsy not considered. 

Following Bejar et al. (2004), I have also studied two 
expressions of exponential decay of the surface density: 
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Finally, I have also investigated the King (1962) profile 
for gravitationally relaxed globular clusters. Close to the 
centre, the surface density can be expressed by: 



cr(r) PS acir) = 



(15) 



H-(r/rc)2' 

where is the core radius and cro is the central surface 
density. In the limit of the cluster, the surface density is: 

(16) 



(j(r-) ^ at(r) = ai ( - - —\ 
\r rtJ 



where rt is the tidal radius (the value of r at which crt(r) 
reaches zero) and cri is a constant. The overall normalised 
cumulative number of objects that embodies crdr) and crtir) 
is, following the nomenclature by King (1962): 

Cartwright & Whitworth (2004) demonstrated that the equal- 
ity is based on a fallacious assumption. The difi'erences between 
(7(r) and 2rp{r), linked to distribution functions of type 2D1 and 
3D2 in the nomenclature by those authors, are, however, too small 
to be considered in this work. 
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Power-law 



where x = {r/rcY, xt = {rt/vc)^ and a::,nax = (rmax/r-c)^. 

Fig. |3] illustrates the fits of /(r) to to evaluate 

the most suitable expression for cr(r). The best match for a 
simple power-law density is aquired for 5 = 0.9. Power-laws 
with 5 S> 1 and 5 <^ 1 provide inaccurate fits. Likewise, 
the exponential profiles cannot predict the large actual sur- 
face density close to the cluster centre (the binned surface 
density profile in Fig. [1] when plotted in logarithmic scale, 
also shows that the innermost bin deviates from the expo- 
nential profile). The overall King profile has the same prob- 
lem. I performed intensive computations, not shown here, 
to cover the {rt,rc) parameter space of the King profile. No 
clear minimum of the exists for the a Orionis radial dis- 
tribution when fitted to the King empirical density law. The 
best solutions were found for all the combinations that sat- 
isfy Tc = 8-12arcmin and S> Tc. The excesses of light at 
large radii of young massive clusters with respect to King 
(and Elson, Fall & Freeman 1987) profile(s), attributed to 
gas expulsion by Goodwin & Bastian (2006) , cannot explain 
the poor fitting for the King profile at small radii in a Ori- 
onis. 

The best general fit is obtained for a composite power- 
law, as shown in Fig. |4] The cumulative number of a Ori- 
onis objects grows proportional to r (i.e. S — 1.0) up to 
^20 arcmin. This size translates into a physical radius of 
~2 pc. At larger separations, /(r) increses with a lower slope, 
indicating an exponent 5 « 0.7. In reality, at such separa- 
tions, an exponential or a limit King [o-t(r)] profile would 
also fit the data. From the extrapolation of the f{r) oc r 
law up to the radius of 30 arcmin, there is deficit of 30-40 
objects in the outermost annulus. As a result, the a Ori- 
onis cluster may be spatially described as a central (r < 
20 arcmin), dense region -the "core"- and an outer (r > 
20 arcmin), more rarified region -the "halo"-. From these 
data, the power-law index transition between the core and 
the halo is quite smooth. However, the relative drop of /(r) 
at r ~ 20 arcmin might be more or less abrupt because of 
our poor knowledge of the a Orionis stellar population at 
large separations from the cluster centre. As an example, 
while more than 90 % of cluster members and candidates 
of the Maynt catalogue in the innermost 10 arcmin have 
known features of extreme youth, this ratio is about 50 % in 
the halo (past spectroscopic, mid-infrared and X-ray anal- 
ysis have been naturally focused on the cluster centre and 
its surroundings). Going deeply into this subject seems to 
be meaningless because of the relatively small amount of 
investigated cluster objects in the halo, contamination by 
young neighbouring stellar populations in Orion (Jeffries 
et al. 2006; Caballero 2007a) or fore-/background stars (e.g. 
Caballero, Burgasser & Klement, in prep.), variable extinc- 
tion to the northeast of the survey area, incompleteness of 
the Mayrit catalogue (Caballero 2007c) and the "anoma- 
lous" radial distribution discussed next. 

The radial distribution of the cluster core is, on the 
contrary to the halo, free of possible systematic errors. The 
/(r) oc r law in the core is translated to surface and volume 
densities o-{r) oc r^^ and p(r) oc r~^. Radial profile investi- 
gations carried out in other very young star-forming regions 
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Figure 3. Normalised cumulative number of a Orionis objects, 
f*{r), in (red) steps plot style and theoretical distribution in 
(blue) solid line. Top panel: power-law distributions for f{r) oc 
j,l/2^ j,2 (^{j-QYa top to bottom); middle panel: exponential 

distributions for a-{r) oc e~'^^ (e = 1/12 arcmin" -'^) and e"""" 
(e = 1/18^ arcmin"^); bottom panel: King distributions for the 
core, limit (dotted) and overall radial profiles (rc = 10 arcmin, rt 
= 200 arcmin). 



Spatial distribution in a Orionis 5 



Power-law in bits 



18M3„|<M<3.7M3<,| 



1 F 



0.4 - 




10 15 20 
r [arcmin] 

Figure 4. Same as Fig. |3] but for a composite theoretical distri- 
bution with /(r) = 1.1 — ) for r < 21 arcsec and /(r) = 

1.0 ^ ^ ) for r > 23 arcsec. 
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Figure 5. Mass-dependent normaUsed cumulative number of 
high-mass cr Orionis stars as a function of the separation to the 
cluster centre. The observed /* (r) for the 340 objects is overplot- 
ted as in Figs. [3] and |4] 



have found similar distributions. Not far away from o Ori- 
onis, Bate, Clarke & McCaughrean (1998) noticed that the 
stars of the Orion Nebula Cluster are distributed with a core 
of uniform volume density and radius rcore = 0.5 arcmin and 
a volume density profile p(r) oc at larger separations. 
Alike power-law volume densitiy distributions have been 
found in other star-forming regions, like Taurus (Fuller & 
Myers 1992; Ward-Thompson et al. 1994), or low-mass cold 
dark molecular clouds, like Barnard 68 (Alves, Lada & Lada 
2001). Furthermore, the power-law index 2 is an "often-used 
initial condition for numerical calculations of star forma- 
tion" (Burkert, Bate & Bodenheimer 1997). The p{r) oc r"^ 
distribution corresponds to a singular, self-gravitating, (ro- 
tating) isothermal sphere. Finally, from fig. 5 in Cartwright 
& Whitworth (2004) and the Q-parameter value for cr Orio- 
nis derived in Section [2. II I estimate that the volume density 
in the cluster varies as p(r) cx j.-i '^±0 '*, xhe -1.7±0.4 in- 
dex is intermediate between those in p Ophiuchus (-1.2) and 
IC 348 (-2.2) and consistent with p(r) oc r~^ . 

The use of the power-law function /(r) oc r, whose cor- 
responding (T(r) diverges at r = 0, has the drawback of 
an incorrect fit in the innermost 1 arcmin of the cluster. 
Characterising the very centre of the cluster is out of the 
scope of this work, since it can be only accomplished with 
high spatial-resolution facilities (e.g. adaptive optics or mid- 
infrared instruments - van Loon & Oliveira 2003; Caballero 
2005). Moreover, the Bate et al. (1998)'s value of rcoro can 
be suitable applied to <j Orionis. There are only seven stars 
at less than 0.5 arcmin to g Ori AB (Caballero 2007b), so 
the central surface density is ctq ~ 9 arcmin"^. This value 
matches well with the actual value of (T(rcore) for /(r) oc r. 
The observational parallelism between the radial distribu- 
tions of the Orion Nebula Cluster and other star-forming 
regions and of cr Orionis evidences that the collapse of an 
isothermal cloud to form a star cluster might be universal. 
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Figure 6. Same as Fig. [S] but for intermediate-mass stars. 

2.3 Mass-dependent radial distribution 

I have also investigated the radial profile of cr Orionis clus- 
ter members and candidates for different mass intervals. 
I have separated the 340 stars and brown dwarfs in four 
mass groups that are equally spaced in logarithmic scale. 
The boundaries between the groups are at about 3.7, 0.77 
and 0.16 Mq. The groups contain, from the most to the least 
massive, 10, 103, 155 and 72 elements. Figs. [5] to [8] show the 
normalised cumulative number of objects for the four groups 
compared to the average distribution. It is manifest that the 
ten stars more massive than 3.7 Mq depart from the general 
trend /(r) oc r, while the brown dwarfs and stars below 
this mass roughly follow it. The massive stars seem to obey 
the King profile close to the cluster centre (with a surface 
density (Tc(r) = ^° — 72-) rather than a power-law, whose 
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Figure 7. Same as Fig. [5] but for low-mass stars. 
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Figure 8. Same as Fig. [S] but for very low-mass stars and 
brown dwarfs. Note the deficit of very low-mass cluster mem- 
ber and candidates at r < 4 arcmin and the inclined raise at r = 
6-7 arcmin 

exponent should be as low as 5 ~ 0.1. Eight of the ten most 
massive a Orionis stars are, besides, within the 10 arcmin- 
radius circle. This agglomeration of early-type stars towards 
the cluster centre is typical in other very young star-forming 
regions, such as the Orion Nebula Cluster, whose centre is 
defined by the OB-type stars of the 9^ Ori multiple sys- 
tem. This result shows again the resemblance between both 
Orion clusters. 

There is no other characteristic feature in the mass- 
dependent distribitution in Figs. [5] to [H] except for a re- 
markable deficit of very low-mass stars and high-mass brown 
dwarfs (0.16 Mq ^ M ^ 0.03 Mq) in the innermost 4 arcmin 
together with a steep raise of /*(r) at 6-7 arcmin (Fig. |8]). 
On the one hand, there are only two representatives of 
this mass interval within 4 arcmin: Mayrit 36273 and the 
Class I object Mayrit 111208. Both of them are the faintest 



sources with /-band photometry identified in the near- 
infrared/optical/X-ray survey in the centre of o Orionis by 
Caballero (2007b). On the other hand, there are almost 20 
low-mass stars and brown dwarfs in the narrow annulus 5- 
8 arcmin. An 0.16 M0-mass object in the cluster, in the up- 
per limit of low-mass interval, has typical magnitudes / ~ 
16.0 mag, J ~ 14.0 mag. These values are far brighter than 
the DENIS and 2MASS completenesses, even in the inner- 
most region affected by the glare of the multiple system 
a Ori. Caballero (2007b) failed to confirm or refute this ab- 
sence of very low-mass stars and high-mass brown dwarfs. 
The raise of /* (r) at 6-7 arcmin indicates, on the contrary, a 
larger density of very low-mass objects in this annulus. One 
can think of many ways which could give rise to the deficit 
of low-mass objects at small radii and excess at intermediate 
radii. Firstly, low-mass objects could actually form in the 6- 
7 arcmin annulus and not in the inner regions (maybe core 
masses were higher in the centre, or competitive accretion 
caused central brown dwarfs to grow into stars). Alterna- 
tively, they could form in the cluster centre, but were ejected 
via dynamical interactions with the massive stars. Low-mass 
objects do not have enough energy, however, to move further 
away from the deep a Ori gravity well of M >5OM0. The 
deficit-excess needs to be explained by theory, but this par- 
ticular set of observations does not give any indication of a 
preferred formation scenario [see Whitworth et al. (2007) for 
a review on the theory of formation of brown dwarfs and very 
low-mass stars]. Further and innovative observations, able to 
avoid the extense glare of the OB system, are required to de- 
termine if the peculiar distribution of very low-mass objects 
in the cluster centre are due to an observational bias or to 
an actual consequence of the formation mechanism. Some 
observational efforts on this topic have been carried out by 
Caballero (2005, 2007b) and Sherry, Walter & Wolk (2005). 

2.4 Azimuthal asymmetry 

Top window in Fig. [9] shows that the distribution of con- 
firmed cluster members and candidates is not whole radi- 
ally symmetric, with an evident lower density to the west 
of a Ori AB with respect to the east. An elongated sub- 
clustering is manifest to the east- northeast of the cluster 
centre, just in the direction to the Horsehead Nebula. On 
the contrary, Bejar et al. (2004) derived that the variation 
of the radial distribution of their very-low mass stars and 
brown dwarfs over their best exponential law fit was Pois- 
sonian, implying no evidence of subclustering. The object 
sample presented in this paper surpasses Bejar et al.'s one 
and allows to corroborate or invalidate their statement. 

I have looked for an azimuthal asymmetry in the a Ori- 
onis cluster in several steps. First, I have generated 1000 
simulated distribution^ following the power law /(r) = 
) . Ten of them (one is highlighted) are shown in 
the bottom window in Fig. The resemblance with the 
actual distribution, in the top window of the Figure, is evi- 
dent. Second, I have divided the survey area in nine regions; 

^ In the general power-law distribution case with index 6, if r is 
a vector of length A^max of uniformly distributed pseudo-random 
radii between r = and Tmax , then the distribution of the vector 
ri/"* follows /(r) = (t;^;^)*- 
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Figure 9. Actual (top) and simulated (bottom) spatial distribu- 
tions of a Orionis objects in the survey area. The coordinate origin 
is in a Ori A. The bottom window shows 10 Monte Carlo sim- 
ulations following /(r) = i^:^ — ) ■ One distribution is marked 
with small filled squares, the remaining ones are marked with 
points. The sizes of the small and big solid circles arc 10 and 
30 arcmin. The eight equal-size segments of annuli used for the 
computation of the asymmetry factor, f , are also indicated. 

a central area, 10-arcmin wide, where the asymmetry is dif- 
ficult to quantify, and eight segments of the annulus in the 
radius interval 10 arcmin < r < 30 arcmin. Third, for each 
simulated distribution, I have computed its corresponding 
asymmetry factor, C, defined by: 

(jVl+jV2)-(jV7+iV8) 

C- , (18) 

^ ' max 

where A'^i and are the numbers of objects of the two most 
populated segments of annulus, and Nj and Ng, axe those of 
the less populated. The mean and the standard deviation 
of the 1000 computed asymmetry factors are C, = 0.0698, 
= 0.0189. The maximum and minimum asymmetry fac- 
tors are (^max = 0.1382 and (^min = 0.0176. The four val- 
ues [C,, (J(^, Cmax and Cmin) are faultless compatible with 
the Poissonian errors within each segment of annulus. Fi- 



nally, I have measured the asymmetry factor for the actual 
o" Orionis distribution: C,* = 0.1853. This value deviates 6.1 
times the from the C, and is significatively larger than 
the i^max among 1000 Monte Carlo simulations. Consider- 
ing as a first order of approximation that the values C^i (i 
= 1...1000) are distributed following a standard normal dis- 
tribution with parameters C, and cr^, then there is a prob- 
ability p = 1 - erf(6.1/2^/^) « 10"^ (where "erf" is the 
error function) that the actual value C," follows such distri- 
bution. Even accounting for generous systematic errors or 
biases, it is highly probable that the radial distribution of 
objects at more than 10 arcmin from the centre of a Orionis 
is azimuthally asymmetric. 

The most populated segments of annulus of the ac- 
tual distribution are, counting anti-clockwise from 12 hours, 
the second (coinciding with the elongated subclustering in 
the direction to the Horsehead Nebula) and the fifth ones. 
Both two segments and the cluster centre spatially coin- 
cide with a filamentary region of maximum emission at the 
12 /xm IRAS passband (see fig. 2 in Oliveira & van Loon 
2004). This "warmer" region has not been identified in works 
based on recent observations with IRAC and MIPS onboard 
Spitzer (e.g. Hernandez et al. 2007) or with the Spatial In- 
frared Imaging Telescope on the Midcourse Space Experi- 
ment satellite (Kraemer et al. 2003). From my data, I cannot 
postulate whether the largest surface density of o Orionis 
objects originally arises from an hypothetical larger density 
of warm dust in the region or, inversely, the filamentary re- 
gion is a consequence of both the low spatial resolution imag- 
ing capabilities of IRAS&nA the largest surface density of ob- 
jects (i.e. the red o Orionis objects, many with mid-infrared 
excesses due to discs -Oliveira et al. 2006; Hernandez et al. 
2007; Caballero et al. 2007-, generate a smooth background 
at 12 /xm that IRAS ^as not able to resolve). 

The accumulation of stars and brown dwarfs in a fil- 
amentary pattern in o Orionis strongly supports some star 
formation scenarios of collapse and fragmentation of a large- 
scale turbulent molecular cloud, especially those that predict 
burst of star formation in filamentary gas (e.g. Bate, Bonnell 
& Bromm 2003) . It is stimulating to notice that these simu- 
lations assumed the contraction of an isothermal, spherical 
molecular cloud with p{r) oc (see Section [2.2p . The fil- 
amentary accumulation in a Orionis is, however, peculiar, 
because no similar arrangements have been found in other 
star-forming regions. For example, Gomez et al. (1993) and 
Larson (1995) found that the subclustering in Taurus- Auriga 
is in the form of star clumps of ~15 components, while Bate 
et al. (1998) showed that in the Orion Nebula Cluster there 
is no subclustering at all. It is obvious that further inves- 
tigations are needed; percolation or two-point correlation 
function of stars are different approaches that can be used. 



3 SUMMARY 

The ~3 Ma-old a Orionis cluster is a perfect laboratory of 
star formation. I have investigated the radial distribution of 
340 cluster members and candidates in a 30 arcmin-radius 
area centred on a Ori AB, taken from Caballero (2007c) . The 
analysis has covered a mass interval from the 18-1-12 M© of 
a Ori AB to the ~0.03 Mq of the faintest brown dwarf de- 
tectable by DENIS. The cluster shows a clear radial density 
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gradient, quantified by the Q-parameter, that accounts for 
the mean separation between members and the EucUdean 
minimum spanning tree of the cluster. I have calculated the 
functional relations between normalised cumulative numbers 
of objects counting from the cluster centre, /(r), and surface 
densities, a{r). Cumulative distribution functions as these 
avoid many problems associated with binning. Among the 
studied radial (power-law, exponential and King) profiles, 
the best fit is for a composite power-law distribution of clus- 
ter members with a core and a rarificd halo. The core extends 
up to ~20 arcmin from the cluster centre and is nicely mod- 
elled by a surface density <T(r) oc r~^, that corresponds to a 
volume density p{r) ot r^^. This volume density matches, in 
its turn, the radial profile in a cluster formed from the col- 
lapse of a self-gravitating, isothermal sphere. The most mas- 
sive a Orionis stars deviate, however, from the general trend 
and are much more concentrated towards the cluster centre. 
There is also an apparent deficit of very low-mass stars and 
high-mass brown dwarfs (0.16 M© > M > 0.035 M©) in 
the innermost 4 arcmin and an excess in the arnmlus at 6- 
7 arcmin to the central Trapezium- like system. Last, there is 
a significant azimuthal asymmetry due to a filament-shape 
overdensity of objects that connects the cluster centre with 
a part of the Horsehead Nebula. This discovery supports the 
formation scenarios that predict burst of star formation in 
filamentary gas. 
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